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Abstract 

We evaluate all two-body decay modes of the heavy scalar tau in the Minimal Su- 
persymmetric Standard Model with complex parameters (cMSSM) and no generation 
mixing. The evaluation is based on a full one-loop calculation of all decay channels, 
also including hard and soft QED radiation. The renormalization of the relevant sec- 
tors is briefly reviewed. The dependence of the heavy scalar tau decay on the relevant 
cMSSM parameters is analyzed numerically, including also the decay to Higgs bosons 
and another scalar lepton or to a tau and the lightest neutralino. We find sizable 
contributions to many partial decay widths and branching ratios. They are mostly of 
0(5 — 10%) of the tree-level results, but can go up to 20%. These contributions are 
potentially important for the correct interpretation of scalar tau decays at the LHC 
and, if kinematically allowed, at the ILC or CLIC. The evaluation of the branching 
ratios of the heavy scalar tau will be implemented into the Fortran code FeynHiggs. 
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1 Introduction 



Beside the Higgs boson search another important task at the LHC is to search for physics ef- 
fects beyond the Standard Model (SM), where the Minimal Supersymmetric Standard Model 
(MSSM) [1] is one of the leading candidates. Two related important tasks are investigating 
the mechanism of electroweak symmetry breaking, as well as the production and measure- 
ment of the properties of Cold Dark Matter (CDM). The most frequently investigated models 
for electroweak symmetry breaking are the Higgs mechanism within the SM and within the 
MSSM. The latter also offers a natural candidate for CDM, the Lightest Supersymmetric 
Particle (LSP), i.e. the lightest neutralino, \i 0- Supersymmetry (SUSY) predicts two 
scalar partners for all SM fermions as well as fermionic partners to all SM bosons. Contrary 
to the case of the SM, in the MSSM two Higgs doublets are required. This results in five 
physical Higgs bosons instead of the single Higgs boson in the SM. These are the light and 
heavy CP-even Higgs bosons, h and H, the CP-odd Higgs boson, A, and the charged Higgs 
bosons, H . In the MSSM with complex parameters (cMSSM) the three neutral Higgs 
bosons mix [3H5], giving rise to the states hi,h2,hs. The tree- level input parameters are the 
charged Higgs boson mass, M#± and tan/3, the ratio of the two vacuum expectation values. 

If SUSY is realized in nature and the scalar quarks and/or the gluino are in the kine- 
matic reach of the LHC, it is expected that these strongly interacting particles are copiously 
produced. The primarily produced strongly interacting particles subsequently decay via cas- 
cades to SM particles and (if i?-parity conservation is assumed, as we do) the LSP. One step 
in these decay chaing^ is often the decay of a scalar tau, fi 2, to a SM particle and the LSP, 
or as a 'competing process' the scalar taus decay to another SUSY particle accompanied 
by a SM particle. Also neutral and charged Higgs bosons can be produced this way. Via 
these decays some characteristics of the LSP and/or Higgs bosons can be assessed, see, e.g., 
Refs. [TIE] and references therein. At any future e + e~ collider (such as ILC or CLIC) a 
precision determination of the properties of the observed particles is expected pHTTj. (For 
combined LHC/ILC analyses and further prospects see Ref. [12] •) Thus, if kinematically 
accessible, the pair production of scalar taus with a subsequent decay to the LSP and/or 
Higgs bosons can yield important information about the lightest neutralino and the Higgs 
sector of the model. 

In order to yield a sufficient accuracy, one-loop corrections to the various scalar tau decay 
modes have to be considered. We take into account all two-body decay modes of the heavy 
scalar tau, fvf, in the MSSM with complex parameters (cMSSM), but we neglect flavor 
violation effects. More specifically, we calculate the full one-loop corrections to the partial 



1 Scalar taus can also be produced directly at the LHC, see for instance Ref. [6], where, however, only 
cross sections for a lighter stau where evaluated numerically. 

2 It should be noted that the purely loop induced decay channels f 1 ~7 have been neglected because 
they yield exactly zero, see Sect. [3] for further details. 
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where x° denotes the neutralinos, the charginos, r and i> r the tau and tau-neutrino 
and Z and W ± the SM gauge bosons. The total decay width is defined as the sum of the 
partial decay widths ([T]) to (JS}, where for a given parameter point several channels may be 
kinematically forbidden. 

As explained above, we are especially interested in the branching ratios (BR) of the 
decays involving a Higgs boson, Eqs. ([!]), (@| as part of an evaluation of a Higgs production 
cross section, or involving the LSP, Eq. fl3]) as part of the measurement of CDM properties 
at the LHC or a future e + e~ collider. Consequently, it is not necessary to investigate three- 
or four-body decay modes. These only play a significant role once the two-body modes are 
kinematically forbidden, and thus the relevant BR's are zero. 

We also concentrate on the decays of f 2 ~ and do not investigate f 2 + decays. In the presence 
of complex phases this would lead to somewhat different results. However, such an analysis 
of CP-violating effects is beyond the scope of this paper. 

Scalar tau decays have been investigated in many analyses over the last decade. Most 
of them were restricted to tree-level evaluations. Existing loop corrections are restricted 
to the MSSM with real parameters (rMSSM). First tree-level results for stau decays in the 
rMSSM were published in Refs. [T3HT5] . Corresponding tree-level results are implemented 
in SDECAY [16]. Tree-level results in the cMSSM can be found in Ref. [T7HT9] . An analysis 
on three-body decays is given in Ref. [2U] . Complete one-loop corrections to sfermion decays 
in the rMSSM involving SM fermions were presented in Ref. [21J. However, no explicit 
numerical results for the full one-loop corrections to stau decays are included in this paper. 
Full one-loop corrections to sfermion decays involving SM gauge bosons in the rMSSM are 
presented in Ref. [22]. However, again no numerical results for scalar tau decays are included. 
One-loop corrections to stau decays in the rMSSM, derived in a pure DR scheme (see below) 
have been made available in the program package SFOLD [23J. 

Several methods have been discussed in the literature to extract the complex parameters 
of the model from experimental measurements. A determination of the trilinear Stau-Higgs 
coupling, A T , in the rMSSM from heavy MSSM Higgs decays was presented in Ref. [21]. 
CP-even observables to extract its phase, ifiA T , have been analyzed in Ref. [2S]. CP-odd 
observables for this determination are investigated in Ref. [2S] (with more details on the 
specific LHC analysis in Ref. [2Z])- Depending on the realized cMSSM parameter space and 
on some further assumptions on the LHC performance, it seems to be possible to obtain 
limits on, e.g., \A T \ and (fA T - 

In this paper we present for the first time a full one-loop calculation for all two-body 
decay channels of the heavier scalar tau in the cMSSM (with no generation mixing), taking 
into account soft and hard QED radiation. In Sect. |2]we briefly review the renormalization of 
all relevant sectors of the cMSSM. Details about the calculation can be found in Sect. |3j and 
the numerical results for all decay channels are presented in Sect. HI The conclusions can be 
found in Sect. [5j The results will be implemented into the Fortran code FeynHiggs [28H3T] . 
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2 The relevant sectors of the complex MSSM 



All the channels dl]) ([6]) are calculated at the one-loop level, including real QED radiation. 
This requires the simultaneous renormalization of several sectors of the cMSSM. In the fol- 
lowing subsections we briefly review these sectors to make this article self-contained. Details 
about the renormalization of most of the sectors can be found in Refs. [32TI35] . 

2.1 The tau lepton/slepton sector of the cMSSM 

For the evaluation of the one-loop contributions to the decay channels in Eqs. (JI} - (jHJ) a 
renormalization of the scalar tau (f) and r-neutrino (v T ) sector is needed (we assume no 
generation mixing). The stau and tau sneutrino mass matrices Mf and Mf r read 

(Ml+ml + Mlc w {P T -Q T sl) m T X* T \ 

T V rn T X T Ml +ml + Mlc w Q T sl) ' {< } 

M 0t = Ml + ll T c w M 2 z (8) 



with 



X T = A T tan/3 . (9) 



Mf L and Mf R are the soft SUSY-breaking mass parameters, where Mf L is equal for all 
members of an SU(2)l doublet. m T and Q T are, respectively, the mass and the charge of 
the corresponding lepton, i^, denotes the isospin of r/u T , and A T is the trilinear soft- 
breaking parameter. Mz and Mw are the masses of the Z and W boson, c w = Mw/Mz, 
and s w = a/1 — c^. Finally we use the short-hand notations c x = cos(x), s x = sin(x). The 
mass matrix Mf can be diagonalized with the help of a unitary transformation Uf , 

D f = U t M fU t = (< U f =gy. (10) 

The mass eigenvalues depend only on \X T \. The scalar tau masses will always be mass 
ordered, i.e. m#, < m; 



rni 

Tl,2 



\ {Ml + Ml) +m 2 T + \llc w M z (11) 



T \\J[MI — M~ R + Mfap(I* - 2Q T sl)Y + 4m2|X r |2 , 
ml = Ml + llc w M\ . (12) 

2.1.1 Renormalization 

The parameter renormalization can be performed as follows, 

M f M f + 5M f , M, T M Pt + 5M St (13) 
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which means that the parameters in the mass matrix M ? are replaced by the renormalized 
parameters and a counterterm. After the expansion contains the counterterm part, 

5M fll = 5M~ L + 2m T 5m T - M 2 z c w Q T 5s 2 w + (J T 3 - Q T s 2 J(c 2 p 8M 2 Z + M% 5c 2p ) , (14) 

5M fl2 — (A* — ix tan (3) 5m T + m T {6A* - [i 5 tan (3 - tan (3 5fi) , (15) 

SM f21 = 5M* fi2 , (16) 

5M f22 = 5M 2 R + 2m T 5m T + M 2 z c 2/3 Q T 5s 2 w + Q T s 2 w (c w 8M 2 + M 2 Z 6c 2p ) , (17) 

5M, r = 5M 2 L + Il r (c 2/3 5M 2 Z + M 2 Z 5c w ) . (18) 

Another possibility for the parameter renormalization of the staus is to start out with 
the physical parameters which corresponds to the replacement: 

U f M f Uj^U f M f Ut + U f 5M f Ut= (< <) + («f /<) (19) 

where Sm^ and Sm 2 2 are the counterterms of the stau mass squares. 5Y T is the countert- 
ernj^l to the stau mixing parameter Y T (which vanishes at tree-level, Y T = 0, and corresponds 
to the off-diagonal entries in D f = XJ f Mf Tjt , Eq. f lTUj) ). Using Eq. f lT9|) one can express 
<5Mf by the counterterms 5m 2 , 5m 2 2 and SY T . Especially for 51SA fl2 one finds 

5M fl2 = U* fll U fl2 (5m 2 1 - 5m 2 2 ) + U? u U f22 5Y T + U fl2 U* f2i 5Y* . (20) 



Eqs. ( floj) and ( 1201) yield a relation between 5Y T , SA T and 5m T , see below. 
For the field renormalization the following procedure is applied, 

g)^ (1+i ,z,)g) with«^(g : , 

^ (1 + \8Z 0r ) ~v T . (22) 

This yields for the renormalized self-energies 

t fll (p 2 ) = Z fll (p 2 ) + \{p 2 - m 2 fl )(6Z fll + 5Z* f J - 5m 2 fi , (23) 

± fl2 (p 2 ) = S fl2 (p 2 ) + \(p 2 - ml)SZ fl2 + \(p 2 - m 2 2 )5Z; 2i - 5Y T , (24) 

t f21 (p 2 ) = S f21 (p 2 ) + \{p 2 - m 2 fi )5Z; i2 + \{p 2 - m 2 2 )5Z f21 - 5Y* , (25) 

t f22 (p 2 ) = S f22 (p 2 ) + \(p 2 - m 2 2 )(5Z f22 + 5Z? 22 ) - 5m 2 2 , (26) 

U(P 2 ) = ^AP 2 ) + UP 2 ~ ml){8Zo r + SZl) - 5m 2 T . (27) 

In order to complete the tau lepton/slepton sector renormalization also for the corresponding 
lepton (i.e. the r mass, m T , and the lepton fields tl, t^, v tl ) renormalization constants have 
to be introduced: 

m T — > m T + 5m T , (28) 



3 The unitary matrix Uf can be expressed by a mixing angle and a corresponding phase. Then the 
counterterm 8Y T can be related to the counterterms of the mixing angle and the phase (see Ref. [36]). 
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Kip) = i ^£ l Ut ( p 2 ) , 




where the components are given by 








sf (p 2 ) = sf (p 2 ) - 


^(8Z^ + 8Zf) - 8m T 




^(8Z R + 5Z^)-5m T 


K T (P 2 ) = ^(P 2 ) + 


\{5Z^+5Z£), 



TL/R ~+ (1 + \&Z L T I R ) T L/R , (29) 

v Tl ^(\ + \5Z Vt )u Tl , (30) 

J I R 

with 5m T being the tau mass counterterm and 5Z T being the Z factors of the left / right- 
handed charged lepton fields; 5Z Ut is the neutrino field renormalization. Then the renormal- 
ized self energy S T can be decomposed into left /right-handed and scalar left /right-handed 
parts, Er and Y,r L ^ SR , respectively, while only the left-handed part exists for the self 
energy £„ r of the massless neutrino 

trip) =iu± L T {p 2 ) (p 2 ) + u;_Sf (p 2 ) + co + ± s T R (p 2 ) , (31) 

(32) 



(33) 
(34) 
(35) 

(36) 

and u± = |(1± 75) are the right- and left-handed projectors, respectively It should be 

noted that Re^ R {p 2 ) = (Re£f L (p 2 ))* holds due to CVT invariance. Re denotes the real 
part with respect to contributions from the loop integral, but leaves the complex couplings 
unaffected. 

2.1.2 The tau neutrino/sneutrino sector 

We follow closely the renormalization presented in Ref. |32l|33], slightly modified to be ap- 
plicable to the tau/stau sector. 

(i) The tau neutrino is defined on-shell (OS), yielding the one- loop field renormalization 

Re5Z UT = -ReS, r (0) , lm8Z Ur = . (37) 

(ii) The v T mass is defined OS, 

ReE, T (m? T ) = 0. (38) 
This yields for the tau sneutrino mass counter terms 

5ml=ReX* T (ml) . (39) 

(iii) Due to m Vr = no off-diagonal parameters in the sneutrino mass matrix have to be 
renormalized. 
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(iv) The diagonal tau sneutrino Z factor is determined OS such that the real part of the 
residuum of the propagator is set to unity, 

ReK T (p 2 )\ p2=ml = . (40) 

with E'Qo 2 ) = 9 - . This condition fixes the real part of the diagonal Z factor to 

Re^ = -Re£',> 2 )| , (41) 

y Or 

which is correct, since the imaginary part of the diagonal Z factor does not contain 
any divergences and can be (implicitly) set to zero, 

WZ* T = 0. (42) 

Including absorptive parts of self-energy type corrections into this Z factor leads to new 
combined factors Z which are (in general) different for incoming particles/outgoing 
antiparticles (unbarred) and outgoing particles/incoming antiparticles (barred), see 
Refs. [33], IM] f° r more details. The combined diagonal tau sneutrino Z factors read 

6Z 0t = -£^> 2 )L 2 , 8Z, t =5Z 9t . (43) 

1 Vr 

(v) Due to m VT e no off-diagonal field renormalization for the tau sneutrino has to be 
performed. 

2.1.3 The tau/stau sector 

We choose the stau masses m^, and the tau mass m T as independent parameters. 
Since we also require an independent renormalization of the scalar neutrino, this requires 
an explicit restoration of the SU(2)l relation, achieved via a shift in the Mf L parameter 
entering the f mass matrix (see also Refs. [371138] )• Requiring the SU(2)l relation to be 
valid at the loop level induces the following shift in M~ (f ) 

M 2 fL (f ) = Ml (y T ) + 5 Ml (y T ) - 5 Ml (f ) (44) 

with 

5Ml(f) = \U fll \ 2 5m\ + \U fl2 \ 2 5ml - U, 22 U;jY T - U Tl2 U;jY T * - 2m T 5m T 

+ Mf c 2f s Q T 6s 2 w - (I 3 T - Q T sl) (c 2/3 5M 2 + Mf «5c 2/3 ) , (45) 

SMl {v T ) = Sm 2 T - II (c 2/3 5M 2 Z + M 2 5c 2/3 ) . (46) 

This choice avoids problems concerning UV- and IR-finiteness as discussed in detail in 
Ref. [52], but also leads to shifts in both stau masses, which are therefore slightly shifted 
away from their on-shell values. An additional shift in M? R recovers at least one on-shell 
stau mass, which is now compatible with our choice of independent parameters 

M 2 (n) = , m 2 |A;-^tan/3| 2 2 _ 2 2 2 

R M? L (l)+ml + Mlc 2 ,(P T -Q T sl)-ml zW^ Tvi n K) 
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The choice of stau for this additional shift, which relates its mass to the stau parameter Mf fl , 
also represents a choice of scenario, with the chosen stau having a dominantly right-handed 
character. A "natural" choice is to preserve the character of the staus in the renormalization 
process. With our choice of mass ordering, m fl < m? 2 (see above), this suggests to recover 
for M~ L > M~ R , and to recover m^ 2 for the other mass hierarchy. Consequently, for our 
numerical choice given below in Tab. CD, we insert m f2 into Eq. (T4"7"|) and recover its original 
value from the re-diagonalization after applying this shift. 

For the tau/stau sector we can now employ a "full" on-shell scheme, where the following 
renormalization conditions are imposed: 

(i) The tau mass is defined on-shell, yielding the one-loop counterterm Sm T : 

Sm T = ±Re {m T [^(m 2 T ) + Z?{m 2 r )] + [£f (m 2 T ) + Zf R {m 2 T )] } , (48) 

referring to the Lorentz decomposition of the self energy E T (p), see Eq. (13 ip . 
The field renormalization constants are given by 

8Z? R = -F^{S^K) + ml (Ef K) + Sf K)) + m T (sf (m 2 T ) + £? R 'K)) 
±^-(SfK)-S T ^K))} (49) 
with S'(m 2 ) = 2 2 . 

\ / Dpi \p2.— m A 

(ii) The stau masses are also determined via on-shell conditions (291139] . yielding 

5ml = neE^«) (< = 1, 2) . (50) 

(iii) The non-diagonal entry of Eq. (fl9l) is fixed as [32"ll3"9| l4T)] 

5Y T = lRe{Z^ 2 (ml) + Z fl Ml)} > (51) 

which corresponds to two separate conditions in the case of a complex 5Y T . The 
counterterm of the trilinear coupling 5A T can be obtained from the relation of Eqs. (fT51) 
and (GOD, 



<L4 T = — [E/^ C/? 2 (5m? - 5m 2 f2 ) + f/ fll U* f22 5Y* + £7? 2 £7 f21 5Y T - (A T - // tan p) 5m T ] 

+ (5 ff tan + (f6tan/3 ) . (52) 

So far undetermined are £tan/3 and 5fi, which are defined via the Higgs sector and the 
chargino/neutralino sector, see Ref. [33] for details. 

(iv) The diagonal scalar tau Z factors are determined OS such that the real parts of the 
residua of the propagators are set to unity, 

Re±' f Jp 2 )\ p2=ml = (2 = 1,2). (53) 
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This condition fixes the real parts of the diagonal Z factors to 

Re SZ fu = -ReZ> fu (p 2 )\ (i = 1,2) , (54) 

which is correct, since the imaginary parts of the diagonal Z factors does not contain 
any divergences and can be (implicitly) set to zero, 

lm5Z fi . = (i = 1,2) . (55) 

Including absorptive parts of self-energy type corrections into these Z factors leads to 
new combined factors Z 

8Z fii = (p 2 ) \ 2 p2=mh , 5Z fii = 5Z fii . (56) 

(v) For the non-diagonal Z factors we impose the condition that for on-shell staus no 
transition from one stau to the other occurs, 

ReE fl2 (m? ) = , ReE f21 (m? ) = (i = 1 , 2) . (57) 

This yields 

ReS fl2 (m? 2 ) -5Y T _ ReEf^mfJ — 5Y* 

"Zfl2 — +2 7 2~^\ ' "'"Til — — 2 5 T\ ' V*°) 

(mi — m% ) (mi — mi ) 

Taking the absorptive parts of the self-energy type corrections into account, the con- 
ditions change to 

E fl2 (m?) = 0, E f21 (m?) = (i = l,2). (59) 

This yields the following combined field renormalization constants for incoming par- 
ticles/outgoing antiparticles (unbarred) and outgoing particles/incoming antiparticles 
(barred), 

E fl2 (m? ) -5Y T - Z f21 (m 2 ~) - 5Y* . 



£f 21 {rn 2 fl ) - 8Y* - E fl2 (m? ) - 5Y T 

(mi — mi ) (m| — mi ^ 



2.2 The Higgs and gauge boson sector of the cMSSM 

The two Higgs doublets of the cMSSM are decomposed in the following way, 

Besides the vacuum expectation values v\ and V2, in Eqs. (162]) a possible new phase £ between 
the two Higgs doublets is introduced. The Higgs potential Vh can be written in powers of 
the Higgs fields, 

V H = . . . + Tfa 0i + T 02 02 + T X1 xi + T X2 \i 
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Hxx 



02 

Xi 

W 



'1 >' 



■)m; 



+ . . . 



(63) 



where the coefficients of the linear terms are called tadpoles and those of the bilinear terms 



are the mass matrices M 



<t><t>xx 



and M„ 



After a rotation to the physical fields one obtains 



V, 



H 



T h h + T H H + T A A 



\ (h, H, A, G) M 



diag 
hHAG 



H 

A 



H~ 
G~ 



(64) 



where the tree-level masses are denoted as m h , m H , rriA, rn G , M H ±, m G ±. With the help 
of a Peccei-Quinn transformation [H] \i and the complex soft SUSY-breaking parameters in 
the Higgs sector can be redefined [12] such that the complex phases vanish at tree-level. 

Concerning the renormalization we follow the usual approach where the gauge-fixing 
term does not receive a net contribution from the renormalization transformations. As 
input parameter we choose the mass of the charged Higgs boson, M#± . All details can be 
found in Refs. [311133^1 (see also Ref. [13] for the alternative effective potential approach and 
Ref. [44] for the renormalization group improved effective potential approach including Higgs 
pole mass effects). 

Including higher-order corrections the three neutral Higgs bosons can mix 



(h, H, A) — > {h x MM) 
where we define the loop corrected masses according to 



M hl < M h2 < M h3 



(65) 



(66) 



A vertex with an external on-shell Higgs boson h n (n = 1,2,3) is obtained from the decay 
widths to the tree-level Higgs bosons via the complex matrix Z [31], 



T/i„ — [ZjniT/i + [Z] n2 IV + [zj^r^ + . . . 



(67) 



where the ellipsis represents contributions from the mixing with the Goldstone boson and 
the Z boson, see Sect. [3j It should be noted that the 'rotation' with Z is not a unitary 
transformation, see Ref. [31] for details. 

Also the charged Higgs boson appearing as an external particle in a stau decay has to 
obey the proper on-shell conditions. This leads to an extra Z factor, 



'H-H+ 



l + Re£' H - H+ (p 2 



2 =M? 



(68) 



4 Corresponding to the convention used in FeynArts/FormCalc, we exchanged in the charged part the 
positive Higgs fields with the negative ones, which is in contrast to Ref. [31]. As we keep the definition of 
the matrix M^,±^± used in |31) the transposed matrix will appear in the expression for M^± g G± . 
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Expanding to one-loop order yields the Z factor that has to be applied to the process with 
an external charged Higgs boson, 



with 



Z H - H+ = 1 + ±6Z H - H+ (69) 



5Z H - H+ = -Re£' H „ H+ (p 2 )\ p2=Ml± = -ReZ' H - H+ (M 2 H± ) - 5Z H - H+ . (70) 



As for the neutral Higgs bosons, there are contributions from the mixing with the Goldstone 
boson and the W boson. This Z factor is by definition UV-finite. However, it contains 
IR-divergences that cancel with the soft photon contributions from the loop diagrams, see 
Sect. EJ 

For the renormalization of tan (3 and the Higgs field renormalization the DR scheme is 
chosen [3TJ[33]. This leads to the introduction of the scale [Ir, which will be fixed later to 
the mass of the decaying particle. 

2.3 The chargino/neutralino sector of the cMSSM 

The mass eigenstates of the charginos can be determined from the matrix 

x= ( M 2 V2sm(3M w \ 

\V2cos/3M w fi J' 1 ' 

In addition to the higgsino mass parameter \i it contains the soft breaking term M 2 , which 
can also be complex in the cMSSM. The rotation to the chargino mass eigenstates is done by 
transforming the original wino and higgsino fields with the help of two unitary 2x2 matrices 
U and V, 



XT=(^) with # = ^QQ. and ^ = ^'U+) ' ( 72) 



where the ith mass eigenstate can be expressed in terms of either the Weyl spinors ipf and 
^ or the Dirac spinor \i ■ These rotations lead to the diagonal mass matrix 

Mj- = V*X T U f = diag(m-±,m-±) . (73) 

From this relation, it becomes clear that the mass ordered chargino masses m-± < m-± 
can be determined as the (real and positive) singular values of X. The singular value 
decomposition of X also yields results for U and V. 

A similar procedure is used for the determination of the neutralino masses and mixing 
matrix, which can both be calculated from the mass matrix 



Y 



/Mi -M z s w cos $ M z s w sin /3 \ 

M 2 M z c w cos -M z c w sin /3 

-M z s w cos (3 M z c w cos —ji 

\ M z s w sin P —M z c w sinP —y, J 



(74) 
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This symmetric matrix contains the additional complex soft-breaking parameter M\. The 
diagonalization of the matrix is achieved by a transformation starting from the original 
bino/wino/higgsino basis, 



xl 



( 



V 




(75) 



(76) 



where denotes the two component Weyl spinor and xl the four component Majorana 
spinor of the fcth neutralino field. The unitary 4x4 matrix N and the physical neutralino 
masses result from a numerical Takagi factorization of Y. The symmetry of Y permits the 
non-trivial condition of using only one matrix N for its diagonalization, in contrast to the 
chargino case shown above. 

Concerning the renormalization we use the results of Ref. [331S5H3S] • Since the chargino 
masses m ~± , m -± and the lightest neutralino mass m f o have been chosen as independent 
parameters the one-loop masses of the heavier neutralinos are obtained from the tree-level 
ones with the shifts 



where the renormalized self energies of the neutralino have been decomposed into their 
left /right-handed and scalar left /right-handed parts as in Eq. ( 13~l~j) . Am^o = is just the 
real part of one of our renormalization conditions. Special care has to be taken in the regions 
of the cMSSM parameter space where the gaugino-higgsino mixing in the chargino sector is 
maximal, i.e. where // ~ Mi- Here 6M2 (see Eq. (180) in Ref. [33]) and S/j, (see Eq. (181) 
in Ref. [33]) diverge as (^11^22^11^2*2 ~~ ^12 ^21 ^12^2*1 ) _1 an d the loop calculation does not 
yield a reliable result. An analysis of various renormalization schemes was recently published 
in Ref. [1H], where this kind of divergences were discussedH In Ref. [1H] it was furthermore 
emphasized that in the case of the renormalization of two chargino and one neutralino mass 
always the most bino-like neutralino has to be renormalized in order to find a numerically 
stable result (see also Ref. [50]). In our numerical set-up, see Sect. HI the lightest neutralino 
is nearly always rather bino-like. If required, however, it would be trivial to change our 
prescription from the lightest neutralino to any other neutralinoj^l 

Similar divergences appearing in the on-shcll renormalization in the sbottom sector, occurring for 
"maximal sbottom mixing", have been observed and discussed in Refs. [32ll33j . 

6 In Ref. [49] it was also suggested that the numerically most stable result is obtained via the renormal- 
ization of one chargino and two neutralinos. However, in our approach, this choice leads to IR divergences, 
since the chargino mass changes (from the tree-level mass to the one-loop pole mass) by a finite shift due to 
the renormalization procedure. Using the shifted mass for the external chargino, but the tree-level mass for 
internal charginos results in IR divergences. On the other hand, in general, inserting the shifted chargino 
mass everywhere yields UV divergences. Consequently, we stick to our choice of imposing on-shell conditions 
for the two charginos and one neutralino. 




(k = 2,3,4) , 



(77) 
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3 Calculation of loop diagrams 



In this section we give some details about the calculation of the higher-order corrections to 
the partial decay widths of scalar taus. Generic diagrams are shown in Figs. [I] - [6j Not 
shown are the diagrams for real (hard or soft) photon radiation. They are obtained from the 
corresponding tree-level diagrams by attaching a photon to the electrically charged particles. 
The internal generically depicted particles in Figs. [I] - E] are labeled as follows: F can be a 
tau r, tau- neutrino z/ r , chargino xf or neutralino x°ki & can be a sfermion fa or a Higgs boson 
h n , V can be a photon 7 or a massive SM gauge boson, Z or W ± . For internally appearing 
Higgs bosons no higher-order corrections to their masses or couplings are taken into account; 
these corrections would correspond to effects beyond one-loop orderly For external Higgs 
bosons, as described in Sect. 12. 2\ the appropriate Z factors are applied and on-shell masses 
(including higher order corrections) are used (as evaluated with FeynHiggs [2"8H3"T] ). 

Also not shown are the diagrams with a gauge-Higgs boson or a Goldstone-Higgs boson 
self-energy contribution on the external Higgs boson leg. They appear in the decay f 2 ~ — > 
f±h n , Fig. [U with a Zj G-h n transition and in the decay fjf —> v T H~ , Fig. [5j with a W~ /G~~- 
H~ transition]^ The corresponding self-energy diagram belonging to the process f 2 ~ — > rf Z 
or f 2 ~ — > i> T W~ , respectively, yields a vanishing contribution for external on-shell gauge 
bosons due to e ■ p = for p 2 = M§ (p 2 = Myy), where p denotes the external momentum 
and e the polarization vector of the gauge boson. 

Furthermore, in general, in Figs. [1] - [6] we have omitted diagrams with self-energy type 
corrections of external (on-shell) particles. While the contributions from the real parts of the 
loop functions are taken into account via the renormalization constants defined by on-shell 
renormalization conditions, the contributions coming from the imaginary part of the loop 
functions can result in an additional (real) correction if multiplied by complex parameters 
(such as A T ). In the analytical and numerical evaluation, these diagrams have been taken 
into account via the prescription outlined in Sect. and their numerical contributions are 
included in the results discussed in Sect. SJ 

Finally it should be noted that the purely loop induced decay channels fjf — > 7 yield 
exactly zero due to the fact that the decay width is proportional to e ■ p and the photon is 
on-shell, i.e. e ■ p — 0. 

The diagrams and corresponding amplitudes have been obtained with FeynArts ver- 
sion 3.7 |51j . The model file, including the MSSM counterterms, is largely based on Ref. |%8] , 
however adjusted to match exactly the renormalization prescription described in Sect. [2j see 
also Refs. [32H35] . The further evaluation has been performed with FormCalc version 7.3 
(and LoopTools version 2.7) [52] . 

Ultraviolet divergences 

As regularization scheme for the UV-divergences we have used constrained differential renor- 
malization [S3], which has been shown to be equivalent to dimensional reduction [SI] at the 
one- loop level [52]. Thus the employed regularization scheme preserves SUSY J55j[56] and 

7 We found that using loop corrected Higgs boson masses in the loops leads to a UV divergent result. 

8 From a technical point of view, the W~ /G~-H~ transitions have been absorbed into the respective 
counterterms, while the Z/G-h n transitions have been calculated explicitly. 
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Ti , - Ti 




Figure 1: Generic Feynman diagrams for the decay — >■ 7"j~/i n (n = 1,2,3). F can be a 
tau, tau-neutrino, chargino or neutralino, S can be a sfermion or a Higgs boson, can be a 
7, Z or W . Not shown are the diagrams with a Z— /i n or G—h n transition contribution on 
the external Higgs boson leg. 



guarantees that the SUSY relations are kept intact, e.g. that the gauge couplings of the 
SM vertices and the Yukawa couplings of the corresponding SUSY vertices also coincide to 
one-loop order in the SUSY limit. Therefore no additional shifts, which might occur when 
using a different regularization scheme, arise. All UV-divergences cancel to all orders in the 
final result. 

Infrared divergences 

The IR-divergences from diagrams with an internal photon have to cancel with the ones from 
the corresponding real soft radiation. They are included via analytical formulas following 
the description given in Ref. [57]. The IR-divergences arising from the diagrams involving 
a 7 are regularized by introducing a photon mass parameter, A. All IR-divergences, i.e. all 
divergences in the limit A — > 0, cancel to all orders once virtual and real diagrams for one 
decay channel are added. 

However, in order to achieve the all order cancellation, special care has to be taken in the 
decay modes involving scalar neutrinos and a W boson. Using tree-level stau masses yields 
a cancellation of IR divergences to all orders for all f% decay modes. However, inserting 
the one-loop corrected stau masses (see Sect. 12. lj ). as required for consistency, we found 
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n n n 





Figure 2: Generic Feynman diagrams for the decay r 2 — >■ f x Z. F can be a tau, tau- 
neutrino, chargino, or neutralino, 5" can be a sfermion or a Higgs boson, V can be a 7, Z or 
W ± . 



cancellation to all orders of the related IR divergences, except for the decay mode f% — > 
v T W~ . Within this decay the tree- level relation required by the SU(2) symmetry M fi (f) = 
Mf L {v T ), corresponding to 

\U fll \ 2 m 2 fl + \Ur 12 \ 2 m 2 f2 = ml T +m 2 T - cos 2(3 , (78) 

has to be fulfilled to yield a cancellation of all IR divergences On the other hand, the 
requirement of on-shell stau masses as well as an intact SU(2) relation at the one-loop level 
leads to the necessity of a shift in the scalar tau masses, see Eq. ( 1441) . Therefore Eq. ( 1781) is 
"violated" at the one- loop level, introducing a two- loop IR divergence in r(f 2 " — > u T W~). 
In order to eliminate this two-loop IR divergence we introduced a counterterm in the f 2 h' T W 
vertex, 

5Z ir = ([|t^| a T4 + \U fl2 \ 2 m 2 2 ] - [\U fll \ 2 ml + |[/ fl2 | 2 r4] shift ) x IR d iv (80) 
9 Eq. (|78| has been deduced via 

M% L (f) = \U fll \ 2 m 2 fl + \U fl2 \ 2 ml_ - M 2 z c 2fS (ll - Q T s 2 w ) - m 2 , (79) 
Ml h {vx)=m\ T -ll T C^M 2 
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Figure 3: Generic Feynman diagrams for the decay fj" — > T~xt (k = 1, 2, 3, 4). F can be a 
tau, tau-neutrino, chargino or neutralino, S can be a stau, tau-sneutrino or a Higgs boson, 
V can be a 7, Z or W ± . 



V T V T V T V T 




Xj Xj Xj 




Xi 



Figure 4: Generic Feynman diagrams for the decay r 2 — > v T Xj (j = 1) 2). F can be a tau, 
tau-neutrino, chargino or neutralino, S can be a stau or a Higgs boson, V can be a 7, Z or 
W ± . 



to restore the tree-level SU(2) relation. The left term in Eq. fl80p contains only 'tree-level' 
values, while the index 'shift' refers to inserting the one-loop masses and mixing matrices. 
The IR divergence has been taken from Eq. (B.5) of Ref. [58] (it can also be found in 
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Figure 5: Generic Feynman diagrams for the decay f% — > i> T H~ . F can be a tau, tau- 
neutrino, chargino or neutralino, S can be a sfermion or a Higgs boson, V can be a 7, Z or 
W ± . Not shown are the diagrams with a W + -H + or G + -H + transition contribution on the 
external Higgs boson leg. 



Ref. ESI), and reads (in our case): 



IR div = * T ln ^ } ^ In ( m \ Mw ) (81) 



27T mf 2 Mvy(l — X 2 ) 

with 



1 - 4m f2 M w /(ml + zO - (M w - m f J 2 ) - 1 

Xt = v T 5 (82) 

*/l - 4m f2 Miy/(m? T + iO - (M w - m f J 2 ) + 1 ' 

where iO denotes an infinitesimally small imaginary part. After including this tree- level 
relation restoring counterterm we find an IR finite results to all orders as required. 

Tree-level formulas 

For completeness we show here also the formulas for the tree-level decay widths: 

I*-(f,f -+ f;K) = |gfff.y.».)l'A"»K,m|„mi.) (n = 3) (g8) 

1 ^^ T ^ Z ^ = 16 „ Ml ml ' (84) 
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r trcc (r^ r-x° fe ) = [ (|C(f 2 -, r-, ^)l| 2 + |C(f 2 ^, r" x° fe )*| 2 ) « - - m|o) 

- 4 Re{ C (f 2 ~ , r ~ , x° ) I C (f 2 ~ , r " , x° ) r } m T m^o 



A 1/2 (m? 2 ,m2,m|o) 



16 7rm| 2 



(& = 1,2,3,4) 



-) = (|C(f 2 -,z/ r ,x7)L| 2 + \C(f 2 -,u T ,xj) R \ 2 ) ( 



~~ m ^±) X 



AV2( m a o,mL ) 



16 7rm| 2 



(i = 1,2) 



r tree (f2- ->• Z> r W") 



16 7r m| 2 



ic^-,^,^-)!^ 3 / 2 ^^^,^ 



16 a M^m 3 f2 



(85) 

(86) 
(87) 
(88) 



where X(x,y,z) = (x — y — z) 2 — 4yz, and the couplings C(a,b,c) can be found in the 
FeynArts model files [60]. C(a, b, c)l,r denote the part of the coupling which is proportional 
to (i=F75)/2. 
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Comparison with other calculations 

As discussed in the introduction, hardly any numerical results for stau decays at the loop 
level are available in the literature. We employed the program SFOLD [23] to obtain nu- 
merical results for scalar tau decays. SFOLD is based on a complete DR renormalization at 
one-loop order (restricted to the rMSSM), but with the possibility also having OS masses 
(instead of the DR masses) internal and/or external^ SFOLD uses a running electromagnetic 
coupling a(Q) with Q denoting the DR scale. This leads to a numerical value significantly 
higher than a(0), see Eq. ( 195]) below. Consequently, our tree- level results differ substan- 
tially. However, at the loop-level the two results are in better agreement as expected. This 
agreement improves for lower values of Q, but differences at the level of 5% were found for 
Q ~ 2 TeV. 

4 Numerical analysis 

In this section we present a numerical analysis of all 12 decay channels. In the various figures 
below we show the partial decay widths and their relative correction at the tree-level ( "tree" ) 
and at the one- loop level ("full"), 

"pfull "ptree 

r tree = r tree (f 2 - xy ) , r fu11 = r fuU (f 2 ~ xy ) , gr/r*™ = rtree , (89) 

where xy denotes the specific final state. The total decay width is defined as the sum of all 
12 partial decay widths, 

pfull "ptree 

cr = E rtree ( f 2~ » r t fu o" = E Tiull ^~ *y) , ^/rsr = • ( 9 °) 

*- tot 

xy xy tot 

We also show the absolute and relative changes of the branching ratios, 



r^fc-^xy) m _ r M (f 2 - -> xy) _ BR fu11 - BR tree 

ftSi ' BK = fMi ' oBR/BR = — M 

1 tot 1 tot i>n 



The last quantity is relevant for an analysis of the impact of the one-loop corrections on the 
phenomenology at the LHC and the ILC 

4.1 Parameter settings 

The renormalization scale [ar has been set to the mass of the decaying particle, i.e. lir = m^ 2 . 
The SM parameters are chosen as follows, see also [6T] : 

• Fermion masses: 

m e = 0.51099891 MeV , m Ue = , 

m M = 105.658367 MeV , m„ M = , 
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It should be noted that we had to use DR masses everywhere for our comparison. 
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m T = 1776.82 MeV , m Vr = , 

m u = 62.8 MeV , m d = 62.8 MeV , 

m c = 1.27 GeV , m s = 101 MeV , 

m t = 172.0 GeV , m b = 4.67 GeV . (92) 

m u and are effective parameters, calculated through the hadronic contributions to: 
A«i? d (Mz) = I £ <#(ln£S-jj)= 0.02793. (93) 

f=u,c,d,s,b f 

• The CKM matrix has been set to unity. 

• Gauge boson masses: 

M z = 91.1876 GeV , M w = 80.399 GeV . (94) 

• Coupling constant: 

a = a(0) = 1/137.035999679 . (95) 

The Higgs sector quantities (masses, mixings, etc.) have been evaluated using FeynHiggs 
(version 2.8.6) [2SM]0 

We will show the results for some representative numerical examples. The parameters are 
chosen according to the scenario S, shown in Tab. [H, but with one of the parameters varied. 
For the scalar quark sector we have chosen M$ L = Mg R = \A q = 1000 GeV (q = u, c, t, d, s, b) 
to yield M hl ~ 120 GeV. The value of Mi is fixed via the GUT relation M x = f tan 2 6» w M 2 w 
|M 2 . The scenarios are defined such that all decay modes are open simultaneously to permit 
an analysis of all channels, i.e. not picking specific parameters for each decay. We will start 
with a variation of m f2 , and show later the results for varying <^>a t - The scenarios are in 
agreement with the MSSM Higgs boson searches at LEP [62],|63], the Tevatron [61] and the 
LHC [65J. Furthermore the following exclusion limits [61] hold in our scenario: 



m^o > 46 GeV, m^o > 62 GeV, m^o > 100 GeV, m^o > 116 GeV, m-± > 94 GeV. (96) 

A few examples of the stau and sneutrino masses in S are shown in Tab. [2j We assume 
SUSY mass scales that allow for the copious production of the colored particles at the LHC, 
with the subsequent cascade decay to uncolored particles we are interested in. Furthermore, 
in S the production of f 2 ~ at the ILC(1000), i.e. with yfs = 1000 GeV, via e + e~ — > f^f^ 
will be possible, with all the subsequent decay modes (JIJ) — dSJ) being open. The clean 
environment of the ILC would permit a detailed study of the scalar tau decays. We find, 
depending on the mixing in the stau sector, cross sections up to a(e + e~ — > fj f 3 _ ) ~ 1 fb, 
where these larger cross sections are found for larger mixing. Even larger cross sections are 
found in the case of a(e + e~ — > f^ff). An integrated luminosity of ~ 1 ab -1 would yield up 
to 1000 scalar taus for a = 1 fb. The ILC environment together with such high numbers of 



As default value within FeynHiggs, /i/? = m t is used. 
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Seen. 


tan (3 


M H ± 


m f2 










4r 


Mi 


M 2 


M 3 


S 


5 


200 


550 


2 m f 2 


1000 


150 


9 


2000 


~ |M 2 


250 


1500 



Table 1: MSSM input parameters for the initial numerical investigation; all masses are 
in GeV. In our analysis Mj l and Mf are chosen such that the values of and rrif 2 are 
realized. For the f sector the shifts in Mf L R (f ) as defined in Eqs. (jUJ) and (I47p are taken 
into account. Mq L R denote the diagonal soft SUSY-breaking parameters in the scalar quark 
mass matrices, while A q is the trilinear squark Higgs coupling, and M 3 denotes the gluino 
mass parameter. The values for Af (/ = T,t,b, . . .) are chosen such that charge- and color- 
breaking minima are avoided [66J. 



Without shifts 


With shifts 


m fl 


rrif 2 










275.000 


550.000 


263.924 


274.478 


550.000 


263.924 



Table 2: The stau and tau-sneutrino masses in S for the numerical investigation; at the 
right-hand side of the table the shifts as defined in Eqs. (144)) and ( 14?]) have been taken into 
account. All masses are in GeV and rounded to one MeV. 



produced staus would result in an accuracy of the relative branching ratio (Eq. fl91j) ) close 
to the statistical uncertainty: a BR of 30% could be determined down to ~ 5%. Depending 
on the combination of allowed decay channels a determination of the branching ratios at the 
few per-cent level might be achievable in the high-luminosity running of the ILC(1000). 

The numerical results we will show in the next subsections are of course dependent on 
choice of the SUSY parameters. Nevertheless, they give an idea of the relevance of the 
full one-loop corrections. As an example, the largest decay widths are r(f 2 ~ — > 7" _ Xi 2 ), 
dominating the total decay width, r tot , and thus the various branching ratios. This is due 
to the strong bino component in Xi 2 m combination with a relatively small mixing in the 
f sector. For other choices of /1, Mi, M 2 , e.g. fi Mi j2 and/or larger mixing in the f sector, 
the light neutralinos would be higgsino dominated and the decay widths would turn out to be 
substantially smaller. Consequently, the corrections to the (other) decay widths would stay 
the same, but the branching ratios would look very different. Channels (and their respective 
one-loop corrections) that may look unobservable due to the smallness of their BR in the 
plots shown below, could become important if other channels are kinematically forbidden. 

4.2 Full one-loop results for varying mf 2 

The results shown in this and the following subsections consist of "tree", which denotes 
the tree-level value and of "full" , which is the partial decay width including all one-loop 
corrections as described in Sect. [3j We start the numerical analysis with partial decay 
widths of f 2 ~ evaluated as a function of m f2 , starting at m f2 = 220 GeV up to m f2 = 2 TeV, 



20 



which roughly coincides with the reach of CLIC. The upper panels contain the results for 
the absolute value of the various partial decay widths, r(f^ — > xy) (left) and the relative 
correction from the full one- loop contributions (right). The lower panels show the same 
results for BR^f^ —> xy). 

Since in this section all parameters are chosen to be real no contributions from absorptive 
parts of self-energy type corrections on external legs can contribute. This will be different 
in Sect. [OI 

In Fig. [7] - M we show the results for the process f 2 — > f^h n {n = 1,2,3) as a function 
of rrif 2 . Here, as well as in the other channels some dips and peaks appear, which are due 
to various thresholds in self-energy or vertex diagram contributions. Three kinks that are 
present in principle in all decays (with the partial exception of f 2 — > u T H~ , see below), 
but that are only partially visible, appear ato m f2 pa 364 GeV ~ Ma + rn fl , 366 GeV ~ 
M H ± +m/> r , 373 GeV ~ Mh + m fl . The thresholds appear in the stau self-energies and thus 
enter via <5m| 2 , [<^Zf] and 5Y T . Visible in the plots is only the kink at m f2 pa 366 GeV. 

In the decays f 2 — > f^h n at nif 2 pa 583 GeV we find a kink due to m fl pa m T + m^o. One 
can see that the size of the corrections of the partial decay widths is especially large very 
close to the production threshold from which on the considered decay mode is kinematically 
possible!^ Away from this threshold relative corrections of ~ +5%, +6%, +6% are found for 
hi,h 2 ,h 3 , respectively. In (all) the plots the value of m f2 for which m fl + m f2 = 1000 GeV 
is shown as a vertical line, i.e. the region where the heavier stau can be produced at the 
ILC(IOOO). In these regions the size of the corrections is only slightly smaller than the 
numbers above. The BRs are at the per-cent level for all three channels. The relative change 
in the BRs for the masses accessible at the ILC(IOOO) are about —5%, —4%, —4% for h±, h 2 , 
h 3 , respectively. For lager masses, only accessible at CLIC, the one-loop corrections are even 
smaller. Depending on the MSSM parameters (and the channels kinematically allowed) the 
one-loop contributions presented here can be relevant for analyses at the ILC and potentially 
as well as at the LHC. 

Next, in Fig. [TDlwe show results for the decay r(f^~ — > rf Z). The dips due to the thresh- 
olds in [SZ?\ 12 22 , SY T and 5m| 2 are the same as before. Furthermore at m f2 pa 416.7 GeV 
a sign change in the stau mixing matrix takes place, and the tree-level result comes out 
zero for Uf = 11. Consequently, around this value the loop corrections are substantially 
larger than the tree- level result, however, not invalidating the perturbative series. The rel- 
ative corrections to the partial decay width in S range between +5% at low m^, i.e. in the 
'TLC(IOOO) regime", to about zero at large m^, with the exception of the region around 
m f2 ^ 416.7 GeV. 

Now we turn to the decays f 2 — > T~Xk = 1> 2, 3, 4), with the results shown in Figs. [11] 
- [Ml Since /i, M\ and M 2 are roughly of the same order, the four states are a mixture of 
gauginos and higgsinos, however, the two lighter states carry a substantial bino component, 
which is the only one in the case of small f mixing, which is not suppressed by small lepton 
masses. Consequently, these two partial decay widths, r(f^" —> T~Xi2) are found to be 
roughly the same and dominating above all other decay widths. Also r(f 2 " — > t~X^a) 

12 Here and below we round most of the values to one GeV. 

13 It should be noted that a calculation very close to threshold requires the inclusion of additional (non- 
rclativistic) contributions, which is beyond the scope of this paper. Consequently, very close to threshold 
our calculation (at the tree- or loop-level) does not provide a very accurate description of the decay width. 
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are roughly the same but significantly smaller, due to their small bino component, where 
the dominating higgsino component is proportional to the (suppressed) Yukawa coupling. 
Apart from the above mentioned general dips and thresholds, we find another threshold 
for r(f^~ —> 1^X2,4) for mf 2 = 300,576 GeV, where m^o = raf 1 + m T , respectively. This 
threshold appears in the stau self-energies and enters via the field renormalization of the 
neutralinos [33]. Small steps can be observed at rrif 2 = 524, 528, 544 GeV in the decay 
f 2 ~ — > t~X4- At this values of the heavy stau mass kinks in the vertex loop function occurs. 
The larger partial decay widths in S for the decay modes f 2 ~ — > T^xt with k = 1,2 go 
up to ~ 5 GeV for rrif 2 = 2 TeV. The radiative corrections are at the 8% (10%) level for 
f 2 ~ — >■ t~x5(2); respectively. Since these decays are dominating the total width the effect in 
the corresponding branching ratios is small. For T(?2 — > t~x^(4)) the corrections are around 
10% (8%), respectively. All four decay modes show a substantial one- loop correction in the 
region accessible at the ILC(IOOO). 

Next in Figs. [151 US] we present the results for f 2 ~ — > v T Xj (j = 1>2). The size of the 
partial decay widths and branching ratios for f 2 ~ ^ u tXi (^2 ~^ v t5G. ) are relatively small, 
below 0.1 (0.01) GeV, respectively. Apart from the general thresholds we find two additional 
ones for each decay, located at m f2 pa 284.8, 287.6 (593.8, 597.0) GeV, where m-± = m T + 

m u T i m x ± = m A* + m ^M ( m x ± = m? ~ m x ± = m > J ' m ' > ^' ^ ne thresholds occur 

in the respective chargino self-energies and thus enter via SM2, 8/j,, 5m^±^ and the field 

renormalization constants. The corrections to the decay widths are ~ +5%, —5% for t 2 ~ — > 
v T X\ and f 2 ~ — > v T X2 m the ILC(1000) relevant region and thus potentially relevant. For 
larger m^ 2 they become smaller in the case of the lighter chargino and larger for the heavier 
chargino. 

We now turn to the decay mode f 2 ~ — > v T H~ , which is shown in Fig. [T71 In addition 
to the general dips and thresholds another one can be found at m f2 pa 601 GeV, where 
mc, T = m T + m-± is realized, i.e. the threshold enters via the sneutrino field renormalization. 
The corrections to r(f 2 ~ — > v T H~) is ~ 4.5% in the ILC(1000) regions, and thus potentially 
relevant, and rises to ~ 8% for large rrir 2 . The relative correction to BR(f 2 _ — > v T H~) reaches 
—5% for low vrif 2 values. 

Finally, results for the other decay mode involving scalar neutrinos, f 2 ~ — > v T W~ , are 
shown in Fig. O We find two additional dips due to thresholds at m f2 = 291, 601 GeV, 
where = m T + m ^f 2 - Furthermore, as in the decay f 2 — > f^Z, at rrif 2 pa 416.7 GeV a 
sign change in the stau mixing matrix takes place, and the tree-level result comes out zero 
for Uf = 1. Consequently, around this value the loop corrections are substantially larger 
than the tree-level result, however, not invalidating the perturbative series (see above). The 
one-loop corrections to the decay width are found to be below the 3% level, except around 
the threshold at pa 416.7 GeV. 
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Figure 7: r(f 2 ~ — > Tree- level and full one- loop corrected partial decay widths 

are shown with the parameters chosen according to S (see Tab. with varied. The 
upper left plot shows the partial decay width, the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR, the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where + = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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Figure 8: r(f 2 ~ — > fih 2 ). Tree- level and full one- loop corrected partial decay widths 
are shown with the parameters chosen according to S (see Tab. with varied. The 
upper left plot shows the partial decay width, the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR, the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where + = 1000 GeV, 
i.e. the maximum reach of the ILC(1000). 
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Figure 9: r(f 2 ~ — > r-f/^). Tree-level and full one-loop corrected partial decay widths 
are shown with the parameters chosen according to S (see Tab. with varied. The 
upper left plot shows the partial decay width, the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR, the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where + = 1000 GeV, 
i.e. the maximum reach of the ILC(1000). 



25 




BR (SBR/BR 




Figure 10: r(f.f — > f^Z). Tree- level and full one- loop corrected partial decay widths 
are shown with the parameters chosen according to S (see Tab. with varied. The 
upper left plot shows the partial decay width, the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR, the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where + = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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Figure 11: r(f^" —> t~Xi)- Tree- level and full one- loop corrected partial decay widths 
are shown with the parameters chosen according to S (see Tab. [Q, with m ?2 varied. The 
upper left plot shows the partial decay width, the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR, the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where + = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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Figure 12: r(f^" —> t^Xq}- Tree-level and full one-loop corrected partial decay widths 
are shown with the parameters chosen according to S (see Tab. [Q, with m ?2 varied. The 
upper left plot shows the partial decay width, the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR, the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where + = 1000 GeV, 
i.e. the maximum reach of the ILC(1000). 
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Figure 13: r(f^" —> t~x®). Tree-level and full one-loop corrected partial decay widths 
are shown with the parameters chosen according to S (see Tab. [Q, with m ?2 varied. The 
upper left plot shows the partial decay width, the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR, the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where + = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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Figure 14: r(f^" —> t~x®). Tree-level and full one-loop corrected partial decay widths 
are shown with the parameters chosen according to S (see Tab. [Q, with m ?2 varied. The 
upper left plot shows the partial decay width, the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR, the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where + = 1000 GeV, 
i.e. the maximum reach of the ILC(1000). 
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Figure 15: r(f^" —> v T Xi)- Tree-level and full one-loop corrected partial decay widths 
are shown with the parameters chosen according to S (see Tab. [Q, with varied. The 
upper left plot shows the partial decay width, the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR, the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where + = 1000 GeV, 
i.e. the maximum reach of the ILC(1000). 
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Figure 16: r(f^" —> v T X2 )• Tree-level and full one-loop corrected partial decay widths 
are shown with the parameters chosen according to S (see Tab. [Q, with varied. The 
upper left plot shows the partial decay width, the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR, the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where + = 1000 GeV, 
i.e. the maximum reach of the ILC(1000). 
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Figure 17: r(f^" —> v T H~). Tree-level and full one-loop corrected partial decay widths 
are shown with the parameters chosen according to S (see Tab. with varied. The 
upper left plot shows the partial decay width, the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR, the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where + = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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Figure 18: r(f^~ — > v T W~). Tree-level and full one-loop corrected partial decay widths 
are shown with the parameters chosen according to S (see Tab. with m ?2 varied. The 
upper left plot shows the partial decay width, the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR, the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where + = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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4.3 Full one-loop results for varying ip^ T 

In this subsection we analyze the various partial decay widths^ and branching ratios as 
a function of (pA T - The other parameters are chosen according to Tab. HJ Thus, within S 
we have + rrif 2 = 825 GeV, i.e. the production channel e + e~ —> f 2 ~ is open at the 
ILC(IOOO). Consequently, the accuracy of the prediction of the various partial decay widths 
and branching ratios should be at the same level (or better) as the anticipated ILC precision. 
It should be noted that already the tree-level prediction depends on <pa t via the stau mixing 
matrix. 

When performing an analysis involving complex parameters it should be noted that the 
results for physical observables are affected only by certain combinations of the complex 
phases of the parameters /x, the trilinear couplings Af (/ = r,t,b,...) and the gaugino 
mass parameters Mi, M 2 , M 3 jl^EZ]. It is possible, for instance, to rotate the phase <pm 2 
away. Experimental constraints on the (combinations of) complex phases arise in particular 
from their contributions to electric dipole moments of the electron and the neutron (see 
Refs. [681169] and references therein), of the deuteron [70] and of heavy quarks [71]. While 
SM contributions enter only at the three-loop level, due to its complex phases the MSSM can 
contribute already at one-loop order. Large phases in the first two generations of sfermions 
can only be accommodated if these generations are assumed to be very heavy [72] or large 
cancellations occur [73], see however the discussion in Ref. [74] . A recent review can be 
found in Ref. [75]. Accordingly (using the convention that </?m 2 = 0, as done in this paper), 
in particular the phase y? M is tightly constrained [76], while the bounds on the phases of the 
third generation trilinear couplings are much weaker. The phases of \x and A T>tt b enter only 
in the combinations (<Pa t t b + ffi) (or in different combinations together with phases of Mi 
or M 3 ). Setting ip^ = (see above) as well as (p^ = (we do not consider this phase in 
this paper) leaves us with the trilinear couplings as the only complex valued parameters. 
The dependence on ipA b and <pA t on the partial decay widths involving scalar bottom and 
top quarks has been analyzed in detail in Refs. [3"2"II3"3~| . and these phases only enter via 
loop corrections into the prediction for the stau decays, whereas A T enters at the tree-level. 
Consequently, we focus on a complex A T and keep A t and Ab real. 

Since now a complex A T can appear in the couplings, contributions from absorptive 
parts of self-energy type corrections on external legs can arise, and they are included in 
the numerical results shown as "full" . The corresponding formulas for an inclusion of these 
absorptive contributions via finite wave function correction factors can be found in Sect. [2J 

As before we start with the decays to Higgs bosons, f% — > rf/i n (n = 1,2,3) shown 
in Fig. [19] - [21] The arrangement of the panels is the same as in the previous subsection. 
In Fig. [191 where the partial decay width r(f 2 ~ — > rf/ii) is given as a function of <pa t , 
one can see that the size of the correction to the partial decay width varies substantially 
with (pA T - The one-loop effects range from +5% to +11% in S. It should be kept in mind 
that the parameters are chosen such that e + e~ — > f± f 2 ~ is kinematically possible at the 
ILC(IOOO) in iS, where the knowledge of such a large variation can be very important. For 
t 2 ~ — > 7"T/i2, shown in Fig. [201 the variation with ipA T is even larger, ranging from +5% 
to +16% with similar conclusions for the ILC(IOOO) as above. The results for f 2 ~ — > f{h^ 

14 Again we note, that we do not investigate the decays of here, which would correspond to an analysis 
of CP-asymmetries, which is beyond the scope of this paper. 
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can be found in Fig. [2TJ Also here the size of the corrections shows a large variation with 
<^>a t i similar to r(f 2 — > ^h 2 ). For the two heavier Higgs bosons, also the variation of 
the respective branching ratios with ipA T is substantial in S, ranging from —4% to +6%, 
potentially exceeding the ILC precision. 

In Fig. [22] we present the phase dependence for the decay mode f 2 ~ — > fyZ. In our 
scenario S the effect of the one-loop corrections to r(f 2 - — > f^Z) varies from ~ —2% to 
~ +10%, again relevant for the ILC precision. An effect of similar size can be observed for 
BR(f 2 _ ffZ). 

In Figs. [231 - [261 we present the variation of r(r 2 ~ — > t~xV) = 1)2,3,4) as a function 
of <pa t - As for the variation with rrif 2 also here for k = 1, 2 larger values of the partial decay 
width are found in S with a similar size as before, again dominating the total width (see the 
discussions above). The one-loop effects on T(f 2 ~ — > t~Xi 2) are about 8%, 10% with a small 
variation with <pa t - F° r r(^T — ^ r ^34); which are substantially smaller, the one- loop effects 
are of similar size, ~ 11%, 7%, respectively, again with a small phase variation. Within S, 
i.e. with the ILC(IOOO) accessible parameter space, the one-loop corrections to the various 
branching ratios is smaller than the effects on the partial widths. However, as discussed 
above, with a different combination of /i, Mi and M 2 the one-loop effects can still exceed 
the potential ILC precision, where, however, only a moderate phase dependence is observed. 

The results for r(f 2 ~ — > v T Xj) (j — 1? 2) are shown in Figs. [271 ES Both decay widths 
change substantially with <pa t - The relative corrections are between +2% and +9% for 
f 2 — )■ u T Xi and between —7% and +3% for f 2 — > v T x^. Within S the variation of the 
branching ratios is slightly smaller for f 2 ~ — > v T Xii and substantially larger, up to —16% for 
r 2 ~ — > v T X2i again with a strong variation with <Pa t , which can then be relevant for the ILC. 

Finally we turn to the decay modes involving scalar neutrinos. In Fig. [221 the results for 
r(f 2 ~ — > v T H~) are presented. The relative correction to the decay width varies strongly 
between +4% and +17% with <pa t , whereas BR^f^ — > v T H~) varies within S between —4% 
and +7%, potentially exceeding the ILC precision. 

The other decay mode involving scalar neutrinos, f 2 ~ — > v T W~ is analyzed in Fig. [301 
The size of the relative correction to the decay width is similar to the f 2 ~ — > v T H~ channel, 
varying between +8% and —2% with <Pa t , with a corresponding variation in the branching 
ratio between and —11%, which is potentially important for physics at the ILC. 
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Figure 19: r(fvT — > f^hi). Tree-level ("tree") and full one-loop ("full") corrected partial 
decay widths (including absorptive self-energy contributions) are shown. The parameters are 
chosen according to S (see Tab. [I]), with <p^ T varied. The upper left plot shows the partial 
decay width, the upper right plot the corresponding relative size of the corrections. The 
lower left plot shows the BR, the lower right plot the relative correction of the BR. 
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Figure 20: r(fvT — > t^/i^)- Tree-level ("tree") and full one-loop ("full") corrected partial 
decay widths (including absorptive self-energy contributions) are shown. The parameters are 
chosen according to S (see Tab. [I]), with <p^ T varied. The upper left plot shows the partial 
decay width, the upper right plot the corresponding relative size of the corrections. The 
lower left plot shows the BR, the lower right plot the relative correction of the BR. 
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Figure 21: r(f^~ — > f^h^). Tree-level ("tree") and full one-loop ("full") corrected partial 
decay widths (including absorptive self-energy contributions) are shown. The parameters are 
chosen according to S (see Tab. [I]), with <p^ T varied. The upper left plot shows the partial 
decay width, the upper right plot the corresponding relative size of the corrections. The 
lower left plot shows the BR, the lower right plot the relative correction of the BR. 
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Figure 22: r(f;f —> f^Z). Tree-level ("tree") and full one-loop ("full") corrected partial 
decay widths (including absorptive self-energy contributions) are shown. The parameters 
are chosen according to S (see Tab. [I]), with <^ T varied. The upper left plot shows the 
partial decay width, the upper right plot the corresponding relative size of the corrections. 
The lower left plot shows the BR, the lower right plot the relative correction of the BR. 
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Figure 23: r(f£" — > r~Xi)- Tree-level ("tree") and full one- loop ("full") corrected partial 
decay widths (including absorptive self-energy contributions) are shown. The parameters are 
chosen according to S (see Tab. [TJ, with <p^ T varied. The upper left plot shows the partial 
decay width, the upper right plot the corresponding relative size of the corrections. The 
lower left plot shows the BR, the lower right plot the relative correction of the BR. 
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Figure 24: r(f£" — > t~j$). Tree-level ("tree") and full one-loop ("full") corrected partial 
decay widths (including absorptive self-energy contributions) are shown. The parameters are 
chosen according to S (see Tab. [TJ, with <p^ T varied. The upper left plot shows the partial 
decay width, the upper right plot the corresponding relative size of the corrections. The 
lower left plot shows the BR, the lower right plot the relative correction of the BR. 
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Figure 25: r(f£" — > T~y$). Tree-level ("tree") and full one- loop ("full") corrected partial 
decay widths (including absorptive self-energy contributions) are shown. The parameters are 
chosen according to S (see Tab. [TJ, with <p^ T varied. The upper left plot shows the partial 
decay width, the upper right plot the corresponding relative size of the corrections. The 
lower left plot shows the BR, the lower right plot the relative correction of the BR. 
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Figure 26: r(f£" — > r^xt)- Tree-level ("tree") and full one-loop ("full") corrected partial 
decay widths (including absorptive self-energy contributions) are shown. The parameters are 
chosen according to S (see Tab. [TJ, with <p^ T varied. The upper left plot shows the partial 
decay width, the upper right plot the corresponding relative size of the corrections. The 
lower left plot shows the BR, the lower right plot the relative correction of the BR. 
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Figure 27: rfr^ —> v T Xi)- Tree-level ("tree") and full one-loop ("full") corrected partial 
decay widths (including absorptive self-energy contributions) are shown. The parameters are 
chosen according to S (see Tab. [T}, with <p^ T varied. The upper left plot shows the partial 
decay width, the upper right plot the corresponding relative size of the corrections. The 
lower left plot shows the BR, the lower right plot the relative correction of the BR. 
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Figure 28: r(f.f —> VrX^)- Tree-level ("tree") and full one-loop ("full") corrected partial 
decay widths (including absorptive self-energy contributions) are shown. The parameters are 
chosen according to S (see Tab. [T}, with <p^ T varied. The upper left plot shows the partial 
decay width, the upper right plot the corresponding relative size of the corrections. The 
lower left plot shows the BR, the lower right plot the relative correction of the BR. 
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Figure 29: r(f.r — > v T H~). Tree-level ("tree") and full one-loop ("full") corrected partial 
decay widths (including absorptive self-energy contributions) are shown. The parameters are 
chosen according to S (see Tab. [T}, with <p^ T varied. The upper left plot shows the partial 
decay width, the upper right plot the corresponding relative size of the corrections. The 
lower left plot shows the BR, the lower right plot the relative correction of the BR. 
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Figure 30: r(f 2 ~ —> v T W~). Tree- level ("tree") and full one-loop ("full") corrected partial 
decay widths (including absorptive self-energy contributions) are shown. The parameters are 
chosen according to S (see Tab. [TJ, with <p^ T varied. The upper left plot shows the partial 
decay width, the upper right plot the corresponding relative size of the corrections. The 
lower left plot shows the BR, the lower right plot the relative correction of the BR. 
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4.4 The total decay width 

Finally we show the results for the total decay width of f 2 _ . In Fig. ED the upper panels 
show the absolute and relative variation with rrtf 2 . The lower panels depict the result for 
varying ipA T - The dips and peaks visible (best) in the upper right panel have been described 
in Sect. 14.21 In S the size of the relative corrections of r tot ranges between about +8% close 
to threshold and goes up to above +9% for large values of m f2 . Such an effect should be 
detectable at the ILC(IOOO) or CLIC. The variation with <£a t is found to be small in our 
numerical scenario, due to the dominance of r(f 2 ~ — > t~x® 2 )i which shows a small variation 
with <Pa t , see Sect. 14.3] The overall size of the effect of (Pa t , shown in the lower row, are 
around +9%, again a value that should be detectable at a future LC. 

5 Conclusions 

We evaluate all partial decay widths corresponding to a two-body decay of the heavy scalar 
tau in the Minimal Supersymmetric Standard Model with complex parameters (cMSSM). 
The decay modes are given in Eqs. (jTJ) — (JHD - The evaluation is based on a complete one- 
loop calculation of all decay channels, also including soft and hard QED radiation. Such 
a calculation is necessary to derive a reliable prediction of any two-body decay branching 
ratio. Three-body decay modes can become sizable only if all the two-body decay channels 
are kinematically (nearly) closed and have thus been neglected throughout the paper. 

We first reviewed the one-loop renormalization procedure of the r/f and v T j v T sector (ac- 
cording to the analyses in Refs. [3"2"ll3~4] ) in the cMSSM, which is relevant for our calculation. 
The details for the Higgs boson and chargino/neutralino sector renormalization can be found 
in Ref. [33]. We have discussed the calculation of the one- loop diagrams, the treatment of 
UV- and IR-divergences that are canceled by the inclusion of soft QED radiation. 

Our calculation set-up can easily be extended to other two-body decay modes in the 
cMSSM. 

For the numerical analysis we have chosen a parameter set that allows simultaneously all 
two-body decay modes, i.e. not to maximize any loop effects. The masses of the scalar taus 
in these scenarios are 275 and 550 GeV for the lighter and the heavier stau, respectively. The 
production of colored particles at the LHC lead to the subsequent cascade decay production 
of scalar taus at the LHC. A decay of the heavy stau to a lighter stau (or sneutrino) and 
a neutral (or charged) Higgs boson can serve as a source of Higgs bosons at the LHC, 
thus a precise knowledge of stau branching ratios is desirable. The scenario also allows 
f ] f f 2 _ production at the ILC(IOOO) or at CLIC, where statistically dominated experimental 
measurements of the heavy stau branching ratios will be possible (depending on the details 
of the MSSM parameters). Depending on the integrated luminosity a precision at the few 
per-cent level could be achievable. 

In our numerical analysis we have shown results for varying rrif 2 and <£a t , the phase 
of the trilinear coupling A T . In the results with varied rrif 2 only the lighter values allow 
fi?2 production at the ILC(IOOO), whereas the results with varied <£a t have sufficiently 
light scalar taus to permit e + e~ — > f^f% . In the numerical scenario we compared the 
tree-level partial widths with the one-loop corrected partial decay widths. In the analysis 
with if>A T varied we explicitly included the effect of the absorptive parts of self-energy type 
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Figure 31: r tot . The tree level ("tree") and full one-loop ("full") corrected total decay 
widths shown with the parameters chosen according to S (see Tab. [1]). The upper left plot 
shows the total decay width, the upper right plot the corresponding relative size of the total 
corrections, with varied. The vertical lines indicate where rrif 2 + = 1000 GeV, i.e. 
the maximum reach of the ILC(1000). The lower plots show the same but with tfA T varied 
(including absorptive self-energy contributions). 
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corrections on external legs. We also analyzed the relative change of the partial decay widths 
to demonstrate the size of the loop corrections on each individual channel. In order to see the 
effect on the experimentally accessible quantities we also show the various branching ratios 
at tree-level (all channels are evaluated at tree-level) and at the one-loop level (with all 
channels evaluated including the full one-loop contributions). Furthermore we presented the 
relative change of the BRs that can directly be compared with the anticipated experimental 
accuracy. 

We found sizable, roughly (9(5 — 10%), corrections in most of the channels. For some 
parts of the parameter space (not only close to thresholds) also larger corrections up to 
15% or even up to 20% have been observed. The size of the full one-loop corrections to the 
partial decay widths and the branching ratios also depends strongly on <$a t - The one- loop 
contributions, again being roughly of (9(5 — 10%), often vary by a factor of 2 as a function 
of ifA T • All results are given in detail in Sects. 14.21 and 14.31 

The numerical results we have shown are, of course, dependent on the choice of the MSSM 
parameters. Nevertheless, they give an idea of the relevance of the full one-loop corrections. 
The largest partial decay widths are r(f 2 ~ — > t~X?2) m our scenario, dominating the total 
decay width, r tot , and thus the various branching ratios. This is due to the strong bino 
component in Xi2- For other choices of /i, M 1; M 2 , e.g. /i -C M 12 , the light neutralinos 
would be higgsino dominated and the decay widths would turn out to be substantially 
smaller. Consequently, corrections to the partial decay widths would stay the same, but the 
branching ratios would look very different. Decay channels (and their respective one- loop 
corrections) that may look unobservable due to the smallness of their BR in our numerical 
examples could become important if other channels are kinematically forbidden. 

Following our analysis it is evident that the full one-loop corrections are mandatory for a 
precise prediction of the various branching ratios. The results for the scalar tau decays will 
be implemented into the Fortran code FeynHiggs. 
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